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ABSTRACT 

We test the relation between relative /-mode frequency variation {5v/v) and 
Lagrangian perturbation in the solar radius {6r/r) obtained by Dziembowski & 
Goode (2004) using several pairs of solar models and show that it doesn't hold 
true for any of the model pairs we have used. We attempt to derive a better 
approximation for the kernel linking the relative frequency changes and the solar 
radius variation in the subsurface layers. 

Subject headings: Sun: oscillations - Sun: helioseismology 

1. Introduction 

The solar /-modes which are essentially surface gravity modes, provide a diagnostic 
of flows and magnetic fields present in the near surface region (Murawski & Roberts 1993; 
Rosenthal & Gough 1994; Rosenthal & Christensen-Dalsgaard 1995; Sofia et al. 2005). The 
/-mode frequencies can also provide an accurate measure of solar radius (Schou et al. 1997; 
Antia 1998). The dependence of solar radius on the 11-year activity cycle is still a matter 
of controversy. Different measurements of solar radius have given conflicting results about 
its temporal variations (e.g., Laclare et al. 1996; Noel, 2004; Kuhn et al. 2004; Chapman 
et al. 2008; Djafer et al. 2008). Using the energy budget of the Sun, it is easy to see 
that the radius variation, if any, would be localized only in the near surface regions which 
are well sampled by /-modes. Hence, we can expect the /-mode frequencies to reflect 
these variations in the solar radius. Dziembowski et al. (2001) and Dziembowski & Goode 
(2004) (henceforth DG04) have obtained a relation between the /-mode frequency variations 
and radius variations in the subsurface layers. This relation was later used by Lefebvre & 
Kosovichev (2005, 2007) to show that helioseismic radius varies in anti-phase with solar 
activity in the outer region of the Sun but there is a change in behavior in deeper layers. 
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However to the best of our knowledge this equation (cf., Eq. (1)) has not been verified. 
In this work, we attempt to test this equation using a few pairs of solar models to check if 
the actual frequency differences match the values computed from Eq. (1). 

The paper is organized as follows: In §2 we test the relation between 5vi/vi and La- 
grangian radius perturbation {6r/r) as derived by DG04 using five different solar models. 
While in §3 we derive a different kernel to describe the above relation and test it with the 
ten model pairs. §4 summarizes the main conclusions from this study. 



DG04 obtained the following relation between the relative /-mode frequency change 
and the Lagrangian radius variation: 



where g is the acceleration due to gravity, // is the mode inertia of /-mode with degree / 
and uJi = 2t[Vi is the angular frequency. If we take two solar models having the same radius 
Rq but differing in structure through differences in input physics, like the equation of state 
or the treatment of convective flux, then in the subsurface layers the radius enclosing the 
same mass will in general not be the same. We can use such models to test Eq. ([T]) by 
comparing the frequency differences to that obtained from the integral on the right hand 
side. In that case 5ui is the frequency difference of the linear adiabatic oscillations of degree / 
and order (/-modes) between the models and 5r is the difference in radius at constant mass 
for the same pair of models. We calculate the eigenfrequencies vi and the eigenfunctions 
for all the /-modes with 40 < / < 1000 for all the models listed in Table 1 using a stellar 
adiabatic pulsation code. Solar f-modes have been recently observed up to a degree of 900 
(Rabello-Soares et al. 2008). In these models the mesh spacing near the surface was reduced 
till there was no difference in the results to ensure that numerical accuracy is sufficient to 
describe the high degree /-modes. All these models have a solar radius of 695780 km and 
use OPAL opacity tables (Iglesias & Rogers 1996). The heavy element abundance Z at the 
surface is 0.0179 in all the models except OPALOWZ for which Z = 0.0127. We have also 
used models with different equations of states (EOS), i.e., the OPAL (Rogers et al. 1996; 
Rogers & Nayfonov 2002), MHD (Dappen et al. 1988; Hummer & Mihalas, 1988; Mihalas et 
al. 1988) and CEFF (Eggleton et al. 1973; Christensen-Dalsgaard & Dappen 1992; Guenther 
et al. 1992). For calculating convective flux, we use either the Mixing Length Theory (MLT) 
or the formulation due to Canuto & Mazzitelli (1991) (CM). We then compare the left hand 
side and the right hand side of Eq. ([1]) for all the ten possible pairs of models. Fig. (lb) 



2. Testing frequency-radius relation with Solar models 
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shows the result for the pair OPALCM - OPALMLT. The models OPALCM and OPALMLT 
are very similar to each other in that the relative frequency difference is of the order of 10~^ 
compared to all other pairs where the same difference is ~ 10~^. Even then it may be noted 
from Fig. (lb) that Eq. ([T]) is a poor approximation. 

The inconsistency in Eq. ([1]) is very prominent in Fig. (2b) where we have compared 
models MHD and OPALOWZ (see Table 1). For this pair 5v/v changes sign near I = 80 
where as the integral on the right hand side of Eq. ([T]) doesn't. This is so because 6r/r 
doesn't change sign in the subsurface layers and hence the integral is always positive (see 



Even though we have shown only two of the ten cases for which we have tested Eq. ([T]) , 
we found similar disagreement for the other eight pairs. Hence it is unlikely, that this 
equation is valid for the Sun and there is considerable doubt about the validity of the results 
obtained by inverting Eq. ([1]) to calculate the solar seismic radius variation with time. 



The equations governing linear adiabatic stellar oscillations are given by (cf., Unno et 



where p, p, c and g are the density, pressure, sound speed and acceleration due to gravity in 
the stellar model. It was precisely this eigenvalue problem which was solved by the adiabatic 
pulsation code in §2. 

The variational formulation of this equation (Chandrasekhar 1964) has been used to find 
relation between frequency variations and structure variations (e.g., Dziembowski et al. 1990; 
Antia & Basu 1994; henceforth AB94). These relations have been obtained by linearizing the 
variational formulation about a reference solar model and are extensively used for structure 
inversions. These relations have also been tested by using a pair of solar models in a similar 
manner to our test in §2. 

Since these relations are obtained by linearizing the differences, they are expected to be 
valid between pairs of models where the differences in structure variables like sound speed 
and density are small. The relative frequency difference, Su/u, between two models can be 
written in terms of four integrals, Ji, I2, 13, 14, as defined in Eq. (8) of AB94. For the sake of 



Fig. (2a)). 



3. Deriving the Kernels for 6r 
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cuVe = V(cVV • e + Vp ■ 
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completeness we provide the expressions of Ii and I3 here. 



= - / p(V ■ O^cV' dr, (3) 

p^'^^v ■ (7) (4) 

It appears that for the /-modes in models hsted in Table 1 the major contribution comes 
from the integral Is involving Sp/p and the second major contribution comes from Ji in- 
volving 5c^ /(?. Fig. (3) shows 5v/v, — (/i + /3)/2ci;^ and —I^/2uj'^ as a function of / for 
model combination OPALCM - OPALMLT. The contributions from I2 and are negligible 
supporting the use of Cowling's approximation by DG04. 

As noted above J3 makes the dominant contribution to 5v/v. Let 5,- denote the variation 
at fixed radius in contrast to 5m which denotes the change at constant mass (the Lagrangian 
variation). The integral J3 includes the derivative of 6rp/ p with respect to radial distance. 
We now attempt to express brpj p in terms of bmT- Using Taylor's theorem we can write: 

braP 5rP , . 

= FT' 

P P Hp 

where Hp is the density scale height. Also conservation of mass gives the following relation: 

p r'^ dr 

Combining Eq. ([5]) with Eq. ([6]) we get: 

^ = _1 A (r^A + (7) 
p r"^ dr " Hp 

Using Eq. ^ and an integration by parts on integral /a in Eq. (j4j) we finally have: 

i^^^-^ + K'^^6mrdr. (8) 

Here K3 = r'^pc^^r^ ■ and the primes denote the derivative with respect to r. 

Transforming the kernel Ki of Sc^/c^ in the integral Ji into a kernel of is more 
algebraically involved. Moreover we make the assumption that = {dp/dp)s = p'/p' 
for the subsurface region. This relation may not be true very near the surface where the 
stratification is not adiabatic. Using hydrostatic balance we can write: 

6rC^ jj d f 5rp\ 6rg 
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Since g oc mjr^ ^ the above equation may be transformed in terms of 5^?" using 

^t9 ^ 47rpr^ bmT 
g m r 

However we neglect the term 6rg/g as it is a small fraction of SmT/r since ATcpr^/m varies 
from 10~^ at O.97i?0 to 10~^ near the surface. Finally a couple of integrations by parts are 
required on Ii due to the presence of the derivative in the first term of Eq. to write it 
completely in terms of 6mr as 

^~^ = -£' {^-^ + ^1 } ^rnr dr, (11) 

where Ki = — p(V • C,)'^c^Hpr'^. A word of caution regarding using the kernel in Ii here is 
necessary. These kernels can have errors very near to the surface, r > O.998-R0 because 
in the outer part of the convection zone the temperature gradient is not close to adiabatic 
value, that is required for the relation = p'/ p' to be valid. We have compared the terms 
-h/2uj^ and -{h + h)/2uj'^ with 5v/v in Fig. (lb) for model pair OPALCM - OPALMLT 
and in Fig. (2b) for the pair MHD - OPALOWZ. It is evident that kernels given in Eq. ([8]) 
and Eq. fllip are in much better agreement with 5v/v than that given by DG04 (our Eq. ([T])). 
From Fig. (2b) it can be seen that 5v/v changes sign at I ~ 80. This is because the slope 
of 5p/ p changes sign around r = 0.987?© (dashed line in Fig. (2a)). There also exists sign 
change in the slope beyond r = 0.998/2©, which is manifested as a sign change in 5v/v at 
/ ~ 450. In contrast 5m^/r (solid line in Fig. (2a)) decreases monotonically and hence the 
sign change in5v/v cannot be captured by the right hand side of Eq. (1). In fact no equation 
of this form with a kernel that doesn't change sign can account for the differences in this 
case. 



4. Conclusions 

We have shown that the relation connecting the variation in relative frequency and the 
Lagrangian radius variation obtained by Dziembowski & Goode (2004) does not hold good 
for any of the ten model combinations considered by us. On the other hand, the kernel 
derived by us gives results in which the discrepancies are much smaller for all the model 
pairs considered in this work. It may be noted that all the models considered by us have the 
same solar radius even though they differ in their subsurface variation in radius at constant 
mass (5mr). In fact Sofia et al. (2005), and Lefebvre & Kosovichev (2005) point out that 
the Solar radius change with solar activity must be non-homologous in the subsurface layers. 
However we would like to point out that in that case the idea of solar radius variation and 
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its observational evidence becomes highly ambiguous. It is not clear if the use of such radius 
variation would give any insight into the structural variations in the solar interior. Moreover 
Eq. ([T]) as well as the relation derived in this work, does not include the effect of magnetic 
field. Thus the effect of magnetic field can be mis-interpreted as due to radius change by 
the use of these equations. Of course, the presence of magnetic field will induce structural 
variations, which may be accounted for by these relations, but in addition there will be 
contribution to frequency variation due to the direct effect of magnetic field, through the 
Lorentz force, which cannot be accounted for by such relations. 

REFERENCES 

Antia, H. M. 1998, A&A, 330, 336 
Antia, H. M., & Basu, S. 1994, A&AS, 107, 421 
Canuto, V. M., & Mazzitelh, I. 1991, ApJ, 370, 295 
Chandrasekhar, S. 1964, ApJ, 139, 664 

Chapman, G. A., Dobias, J. J., & Wahon, S. R. 2008, ApJ, 681, 1698 

Christensen-Dalsgaard, J., & Dappen, W. 1992, A&A Rev, 4, 267 

Dappen, W., Mihalas, D., Hummer, D. C, & Mihalas, B. W. 1988, ApJ, 332, 261 

Djafer, D., Thuillier, G., & Sofia, S. 2008, ApJ, 676, 651 

Dziembowski, W. A., & Goode, R R. 2004, ApJ, 600, 464 

Dziembowski, W. A., Goode, R R., & Schou, J. 2001, ApJ, 553, 897 

Dziembowski, W. A., Pamyatnykh, A. A., & Sienkiewicz, R. 1990, MNRAS, 244, 542 

Eggleton, R R, Faulkner, J., & Flannery B. R 1973, A& A, 23, 325 

Guenther, D. B., Demarque, R, Kim, Y.-C, & Pinsonneault, M. H. 1992, ApJ, 387, 372 

Hummer, D. G., Mihalas, D. 1988, ApJ, 331, 794 

Iglesias, C. A. & Rogers, F. J. 1996, ApJ, 464, 943 

Kuhn, J. R., Bush, R. I., Emilio, M., & Scherrer, R H. 2004, ApJ, 613, 1241 
Laclare, F., Delmas, C, Coin, J. R, & Irbah, A. 1996, Sol. Phys., 166, 211 



- 7- 

Lefebvre, S., & Kosovichev, A. G. 2007, ApJ, 658, L135 
Lefebvre, S., & Kosovichev, A. G. 2005, ApJ, 633, L149 
Mihalas, D., Dappen, W., & Hummer, D. G. 1988, ApJ, 331, 815 
Murawski, K., & Roberts, B. 1993, A& A, 272, 595 
Noel, F., 2004, A & A, 413, 725 

Rabello-Soares, M. C., Korzennik, S. G., & Schou, J. 2008, Sol. Phys., 251, 197 

Rogers, F. J., & Nayfonov, A. 2002, ApJ, 576, 1064 

Rogers, F. J., Swenson, F. J., & Iglesias, C. A. 1996, ApJ, 456, 902 

Rosenthal, C. S., & Christensen-Dalsgaard, J. 1995, MNRAS 276, 1003 

Rosenthal, C. S., & Gough, D. O. 1994, ApJ 423, 488 

Schou, J. 1999, ApJ, 523, L181 

Schou, J., Kosovichev, A. G., Goode, R R., & Dziembowski, W. A. 1997, ApJ, 489, L197 

Sofia, S., Basu, S., Demarque, R, Li, L., & Thuillier, G. 2005, ApJ, 632, L147 

Unno, W., Osaki, Y., Ando, H., & Shibahashi, H., 1989, Nonradial Oscillations of Stars, 2nd 
Edition (University of Tokyo Press, Tokyo) 



This preprint was prepared with the A AS IM^jX macros v5.2. 



-8- 



Table 1: Overview of the solar models used 





Model 


EOS 


Convection 


Z 


1 


OPALCM 


OPAL 


CM 


0.0179 


2 


OPALMLT 


OPAL 


MLT 


0.0179 


3 


OPALOWZ 


OPAL 


CM 


0.0127 


4 


MHD 


MHD 


CM 


0.0179 


5 


CEFF 


CEFF 


CM 


0.0179 
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Fig. 1. — Comparison between models OPALCM and OPALMLT. (a) 5r/r (solid line) at 
constant mass as a function of x = t/Rq for the same pair. The dashed line gives 5p/p 
at constant radius as a function of x. (b) (thick solid), right hand side of Eq. (1) 

(dashed dotted), — (Ji + l-^jloj^ (dashed), —I^/2uj'^ (thin solid) as a function of I {see text 
for definitions). The error bars denote the errors for observed /-mode frequencies (Schou 
1999). Note that the result deteriorates by adding Ji to J3. As mentioned in §3, Ji is not a 
good approximation to Ji in the outermost layers where this pair of models have maximum 
differences. 
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Fig. 2.— Same as Fig. 1 but for models MHD and OPALOWZ 
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Fig. 3.— 5u/u (solid) for OPALCM - OPALMLT, summation of integrals (from AB94) 
-(/i + h)/2uo'^ (dashed dotted), -h/2uo^ (dashed) function of /. Compare this with 
right hand side of Eq. (1) (thick dashed) as a function of /. 



